function m2 = calc_m2(spectrum)
K1 = spectrum.K1;
K2 = spectrum.K2;
K3 = spectrum.K3;
K4 = spectrum.K4;
z1 = spectrum.z1;
z2 = spectrum.z2;
z3 = spectrum.z3;
z4 = spectrum.z4;        
m2 = sum(sum(K1.*z1.*z1+K2.*z2.*z2+K3.*z3.*z3+K4.*z4.*z4));
